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ABSTRACT 



^ I Weak lensing by large scale structure induces correlated ellipticities in the 

I images of distant galaxies. The two-point correlation is determined by the 

^ ' matter power spectrum along the line of sight. We use the fully nonlinear 

evolution of the power spectrum to compute the predicted ellipticity correlation. 
We present results for different measures of the second moment for angular 
^ I scales ^ ~ 1' — 3°, and for alternative normalizations of the power spectrum, in 

Q-i| order to explore the best strategy for constraining the cosmological parameters. 

O I Normalizing to observed cluster abundance the rms amplitude of ellipticity 

\ within a 15' radius is ~ 0.01 2°'^, almost independent of the cosmological model, 

^ ' with Zg being the median redshift of background galaxies. 

Nonlinear effects in the evolution of the power spectrum significantly enhance 
^ ' the ellipticity for 6' < 10' — for 6* ~ 1' the rms ellipticity is ^ 0.05, which is 

I nearly twice as large as the linear prediction. This enhancement means that 

the signal to noise for the ellipticity is only weakly increasing with angle for 
2' < 6 < 2°, unlike the expectation from linear theory that signal to noise is 
strongly peaked on degree scales. The scaling with cosmological parameters 
also changes due to nonlinear effects. By measuring the correlations on small 
(nonlinear) and large (linear) angular scales, different cosmological parameters 
can be independently constrained to obtain a model independent estimate of 
both power spectrum amplitude and matter density Q^- Nonlinear effects 
also modify the probability distribution of the ellipticity. Using second order 
perturbation theory we find that over most of the range of interest there are 
significant deviations from a normal distribution. 
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1. Introduction 



Mapping the large-scale structure (LSS) of the universe is one of the major goals of 
observational cosmology. Traditionally this is performed using large surveys of galaxies, 
either the projected 2-dimensional distributions or the 3-dimensional surveys in redshift 
space. The main shortcoming of the galaxy surveys is that they trace light, while most of 
the matter appears to be dark. One therefore needs to translate the galaxy power spectrum 
into the matter power spectrum and in order to do so one has to make some assumptions 
on the nature of galaxy biasing. Even in the simplest model this can be achieved only up 
to an unknown biasing parameter b, which at present cannot be theoretically estimated. 
While sophisticated N-body and hydro-dynamical simulations will eventually provide some 
answers to this question, at present the biasing relation between the light and matter 
remains rather poorly understood and prevents one from drawing definitive conclusions on 
the amplitude and distribution of mass fluctuations in the universe from galaxy survey data. 

It is clearly important to seek ways to estimate the large-scale structure that are 
insensitive to biasing. Several observable tracers have been proposed that probe directly 
the underlying mass distribution: cosmic microwave background (CMB) anisotropies, 
gravitational lensing, peculiar velocity flows and abundances of massive nonlinear objects. 
This paper focuses on tracing the dark matter with gravitational lensing, in particular 
the effects of weak lensing by large scale structure on background galaxies. Weak lensing 
magnifies and shears the images of distant galaxies. The shear induces an ellipticity in the 
image of an intrinsically circular galaxy. Background galaxies are of course not circular, 
but by averaging over the observed ellipticities of a large number of galaxies, the induced 
ellipticity can be measured, and related to the mass fiuctuations along the line of sight 
and to the spatial geometry of the universe. Ellipticities of distant background galaxies 
averaged over several arcminute windows are sensitive to the mass power spectrum on 
scales of 1 — 10h~^ Mpc. For a given spectrum of mass fiuctuations, it is sensitive to the 
cosmological parameters Qm and f^A- Thus while strong lensing which leads to multiple 
images probes non-typical regions of the universe which contain massive halos, weak lensing 
provides a different and a more direct measure of the mass fluctuations on large scales. 

The flrst calculations of the shear signal due to weak lensing that used modern 
models for the large scale structure power spectrum were those of [Blandford et al.| 



1(1991)1 , pVliralda-Escude (1991)| and [Kaiser (1992)| , based on the pioneering work by Gunn 



(1967). Our work generalizes the results of these authors to include the effects of nonlinear 
evolution of the matter fluctuations for flat as well as open and A— dominated cosmologies. 
Villumsen (1996) has considered some aspects of the linear calculation for open models. 
Very recently, Bernardeau et al. (1996), Kaiser (1996) and Stebbins (1996) have also made 
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the linear calculation for different cosmologies. Our work extends the results of the above 
authors by using the nonlinear power spectrum, which makes a significant difference on 
small angular scales of ^ < 10'. For the nonlinear calculation we have used the prescription 
of Hamilton et al. (1991) as implemented by Jain, Mo fc White (1995"] and Peacock 



Dodds (1996)1 which provide accurate power spectra for different models, valid from the 



linear to the strongly nonlinear regime. 

One can take two different approaches in interpreting a possible measurement of the 
shear signal. The first is to work within the framework of a physical model for the dark 
matter and background cosmology. Such a model is best normalized to COBE and then its 
small scale predictions can be compared to the observational constraints of e.g. cluster or 
damped Ly-a system abundances, strong lensing statistics or peculiar velocity flows, or in 
our case the shear amplitude on a given angular scale. Another possibility, which is less 
model dependent, is to compare the constraints from different tracers on the same physical 
scale. This way one can test the gravitational instability assumption that both tracers probe 
the same underlying power spectrum and place constraints on the cosmological parameters 
which scale differently with the tracers. At the present the CMB data constrain the power 
spectrum only on very large scales, where there are no available data from other tracers. 
While several current CMB experiments are approaching the scales probed by other tracers, 
cosmological parameters such as baryon density, reionization epoch and Hubble constant 
become important and complicate the power spectrum reconstruction. The test mentioned 
above is possible between the cluster abundances (or strong lensing statistics of multiple 
images at large separation, both of which trace essentially the same property) and the 
peculiar velocities. Unfortunately this method cannot give model independent constraints 
on cosmological parameters, because the two tracers scale roughly equally with the matter 
density (Peebles 1980, White, Efstathiou & Frenk 1993). 

A similar comparison can also be made between weak lensing observations and the 
tracers discussed above. For example, the shear amplitude on 15' angular scale and rich 
cluster abundance both probe power spectrum scales around Sh~^ Mpc, so the ratio of 
the two is roughly independent of the shape of the power spectrum. One could hope that 
using such comparisons would lead to an estimate of the mean density of the universe, 
because the two tracers scale differently with it. Unfortunately as shown in this paper 
this is not the case so that this is not a promising method of obtaining the mean density. 
Nevertheless, weak lensing can also give important constraints on the power spectrum both 
on smaller and larger scales, where other tracers give less stringent constraints. On the 
smaller angular scales, our calculations provide the enhancement in the amplitude and 
the change in the scalings with as and Qm that arise due to nonlinear evolution. Once 
reliable detections of the shear at different angular scales are available, these results can be 
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used to constrain cosmological parameters and the matter power spectrum. The principal 
advantage of working with weak lensing is that no assumptions about the formation of 
observable structures need to be made to connect the mass power spectrum to the observed 
shear. The only unknown parameter is the redshift distribution of source galaxies, and in 
principle this can also be determined observationally. 

In §2 and the Appendix, the formalism for the weak lensing calculation is presented 
following Seljak (1995, 1996), who generalized the work by [Blandford et al. (1991)| , [Miralda-] 
Escude (1991)1 and [Kaiser (1992)| to a non-fiat universe and non-linear regime. The 
derivation presented here is complementary to recent derivations by Bernardeau et al. 
(1996) and Kaiser (1996). The linear and nonlinear evolution of the power spectrum is 
introduced in §2.1, and the resulting dependence on cosmological parameters for power law 
spectra is obtained in §2.2. In §3 realistic CDM-like spectra are used to predict the rms 
shear. We provide accurate power law fits for the dependence of the shear on the parameters 
6, Zs, cs and VLm for different cosmological models. §3.1 provides predictions of the shear 
for COBE and cluster-abundance normalized spectra for the range of angular scales which 
are being probed by current and forthcoming observations. The effects of non-gaussianity 
in the distribution of shear are considered in §4 where we compute the skewness of the 
distribution. This is followed by a discussion and conclusions in §5. 



2. Theory of weak lensing 



Gravitational lensing shears and magnifies the images of distant galaxies. The relation 
of the shear to perturbations in the gravitational potential along the line of sight is 
developed in the appendix. Here we shall use the formalism presented in the appendix to 
derive expressions for 3 different measures of the second moment of the shear. 

The observable mean ellipticity of galaxy images can be simply defined in the 
approximation that the images are ellipses with complex eccentricity e given by the axes 
lengths h and a as: e = (6^ — a'^)/{b'^ + a?) e^*^, where ip is the position angle of the 
major axis. In the limit of weak distortions, the eccentricity is the same as ellipticity 
e ^ (1 — h/aje^"^^. While individual galaxies have intrinsic ellipticities which cannot be 
separated from gravitational stretching, we have assumed here that averaging over several 
galaxies leaves only the gravitational component which can in principle be measured. This 



mean ellipticity (e) will be denoted as p and is defined by equation (|A6|) in the appendix. 
We shall use p to denote its amplitude in what follows. 
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Information on the ellipticity is contained in the trace-free part of the shear tensor, 
which is obtained by integrating over all the deflectors between us and the galaxy and over 
the distribution of background galaxies W{x) (see appendix), 



9ix) = r{x) r • (1) 



X 



Here is the gravitational potential and x is the radial comoving distance with xo being 
the horizon distance. We also introduced the comoving angular distance r{x) defined in 
equation (|A2|) . The radial distribution of background galaxies is described with the window 
W{x')- If 3^11 the galaxies are assumed to lie at the same redshift Zs, corresponding to xs, 
then W{x') = ^d{x' ~ Xs): and g{x) takes the simple form, 

«W-'-(x)%^. (2) 

There are several two-point statistics that can quantify the induced ellipticites in galaxy 
images. The easiest to obtain from observations is the rms ellipticity p{6) = {p{d)p{9)*Y^'^ , 
where the overbar indicates the average within a circular aperture of radius 9. Equation 
(H) can be used to express p{9) in terms of the power spectrum of density perturbations 
Ps{k,x)- Following the derivation presented in the appendix we obtain for the ellipticity 
variance 

p\9) = 9{2nf QlH',9-' ^^^dx J Ps{l/r9,x)W^mi 

roc fxa 

= 367r2 / kdk / a~\x) P5{k,x) 9'{x)W^[kr{x)9]dx, (3) 
Jo Jo 

where W2{x) = 2Ji{x)/x with Ji{x) being the Bessel function of first order and a is the 
expansion factor normalized to unity today. A related quantity is the two-point polarization 
correlation function Cpp{9) = (e(0) e*{9)), which is given by a similar expression to the one 
above, replacing the term W2 with Jq: 

Cpp{9) = 36n' nl / kdk / a-\x) Psik, x) g'ix) Mkr{x)9] dx- (4) 
JO Jo 

The Fourier transform of Cpp{9) is the angular ellipticity power spectrum P{1), which is 
often the optimal statistic to use (Kaiser 1992). Here we use the quantity cr^(0 = 2ttI'^P{1), 
which gives the contribution to the variance per log-interval in / and is given by 

a'il) = 36vr^ nl f f ^ a-\x) Ps U = x) dx- (5) 
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This gives an estimate of correlations at an angle 9 = l/l. 

No assumption on the matter power spectrum has been made in the above expressions 
and one can use them both in the linear and nonlinear regime. The nonlinear spectrum is in 
general a non-separable function of wavenumber and redshift, thus leading to a complicated 
coupling of the dependences on the distance factors and the growth of perturbations. It also 
leads to a nonlinear dependence of the predicted ellipticity correlations on the shape and 
amplitude of the initial power spectrum. Thus the amplitude of the ellipticity correlations 
as well as their dependence on cosmological parameters can change in the nonlinear regime 
(i.e., on small angular scales). In the rest of the paper we shall explore in detail the outcome 
of these effects of nonlinear evolution. 



2.1. Evolution of the power spectrum 



If the gravitational potential changes in time its power spectrum will depend on the 
radial distance x- linear theory this dependence is independent of k and can be written 
in terms of the potential growth factor F{x) as P(p{k,x) = F"^ {x) P<t>{^) ■> which gives for the 
density power spectrum Pi{k^x) / '^'^ix) = F^{x)Ps{k) = [{D^/a)^]{x)P5{k), where -D+(x) is 
the linear growth factor for the density. This can be approximated as (Lahav et al. 1991) 

Fix) = 2.5 a-^ (xf + 1.5fi„a-^ + ^k)'^ 
x = l + nU^-' - 1) + fiA(a2 - 1) ; / = f . (6) 

We have ignored the weak dependence of the logarithmic growth factor. For a = 1, this 
expression simplifies to 

Fix = 0) = 2.5Qm (1 + ^"J + 0.5 - ^a)"' • (7) 

For a flat = 1 model the gravitational potential does not change in time in the linear 
regime and the region where gix) peaks dominates the radial integral in equations (^^. 
Typically this is at half the mean comoving distance to the galaxies, so if background 
galaxies lie at z = 1 typical deflectors lie at z ^ 0.3. In low models the gravitational 
potential increases with x, differently for open and for cosmological constant models. In the 
small z limit the growth factor only depends on fi^, (Villumsen 1996). 

On scales where A^(A;,x) = AiTk^Ps{k,x) approaches or exceeds unity nonlinear 
evolution of the power spectrum becomes important. This is particularly important if 
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one is discussing low fi^ models, where both cluster abundances and peculiar velocity 
normalizations give higher density normalization on scales of interest. In the quasilinear 
regime the nonlinear spectrum can be computed using perturbation theory. An alternative 
semi-analytic approach which is accurate from the linear to the strongly nonlinear regime 
(up to density contrasts ~ 10^) was proposed by Hamilton et al. (1991) and further 
developed in subsequent work (e.g., Nityananda & Padmanabhan 1994; Peacock & Dodds 
1996; Jain, Mo & White 1995; Padmanabhan et al. 1996). This prescription involves 
mapping the nonlinear spectrum at a wavenumber k to the linear spectrum at a unique 
wavenumber < k which is given by a spherical collapse model. This leads to a mapping 
of the form A'^{k,x) = G [A|(A;l,x)], where the dimensionless power is defined above, 
and G is a function which varies in a simple way for different initial spectra. The linear 
wavenumber k^ is given hj k^ = k [1 + A^(A;, x)]~^^^- The functional form of the mapping 
is thus specified by the function G, which has been calibrated using high resolution N-body 
simulations and is accurate for a wide range of initial spectra. 



For the predictions of the rms ellipticity in §3 we shall use the fitting formulae of Jain 



Mo fc White (1995)| for Q^n = 1 and of [Peacock fc Dodds (1996)| for the open and A— models 



to describe the nonlinear power spectrum as a function of wavenumber and redshift. We 
refer the reader to the above references for details of the formulae and their implementation. 
While the formulae can differ by a few tens of percent for certain spectra, we find that 
our results for the ellipticity are not affected by more than a few percent as the ellipticity 
integrates the power spectrum over a range of wavenumbers and redshifts. In particular the 
results for the CDM-like models shown in Figures 4-9 would be completely unaffected by 
our choice of fitting formulae. 



2.2. Dependence on Zg, erg and fi^ for power law spectra 



In the next section we shall use the nonlinear evolution of CDM-like spectra to predict 
the rms ellipticity signal. Here we consider the simplified case of power law initial spectra 
Ps{k) = Ak"- to compute the scalings of Cpp{6) and p^iO) with source redshift Zg, and the 
cosmo logical model parameters as and Vim- The parameter is the rms fluctuation in the 
mass on scales of 8/;.~^Mpc, and serves to normalize the power spectrum. We use the linear 
evolution of the power spectrum for which some of the scalings can be obtained analytically. 
We also consider the qualitative modification due to nonlinear evolution by using the stable 
clustering regime which gives the maximal effect of nonlinearities. The length scales that 
contribute to the ellipticity are of order 8/i~^ Mpc for 9 ~ 15'. On these scales most realistic 
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spectra have a spectral index n ~ — 1. On smaller angular scales, the contribution to the 
ellipticity comes from higher wavenumbers which sample a steeper part of the spectrum 
(lower n). We therefore choose power law spectra with n = — 1 and n = —2 for comparison 
with realistic spectra on angular scales of interest. 

For a power law spectrum the A;— integral for Cpp is analytic, and gives: 

Cppi9)=36n^nlA X^vi-nm ^ ^^^'^^^ ^xY'^F^x). (8) 

valid for —2 < n < —1/2. The constant A oc (t| gives the normalization of the power 
spectrum. Aside from the numerical prefactors, the result for p'^{9) is the same and involves 
the same x integral. 

For simplicity we shall take all the source galaxies to be at the same redshift - this 
is a good approximation as the result is insensitive to the source distribution for a given 
median source redshift. The window function g{x) then takes the form given in equation 
(0). For VLrn = 1, r{x) = X) ci-nd the % integral is analytic as well giving the following power 
law dependence: 

, (X al x^-" r . (9) 

Using X = 2-^0"^ (1 — o^^^), the dependence of Cpp on the source redshift Zg can be obtained 
from the above equation. 

For Qrn < 1 and > 0, the result for Cpp is not analytic except in the limit of z -C 1 
which is not of interest for realistic situations. It is simple to understand qualitatively the 
two physical effects that enter, {i) For a given redshift, the distance x increases as flm 
decreases and increases. Therefore there is more path length in the line of sight integral 
as xs is larger for a given Zg. The factor g{x) involves factors of the comoving angular 
diameter distance r{x),''^{xs) and r{xs — x)i it also increases with decreasing flm and 
increasing Q\. (ii) The second factor is the linear growth factor F{x)- Growth of structure 
slows down at low redshift in an open universe, and to a lesser extent in a A— dominated 
universe. Therefore, at a given redshift the linear growth factor, normalized to unity today, 
increases as flm decreases, and for a given Qm it decreases as increases. 

The effect of both the distance and growth factors is to increase the contribution from 
the X integral in equation (H) relative to the Einstein-de Sitter case. The net contribution 
increases as flm decreases, and for a given flrn, as fl\ increases. In the case of A, the 
distance factors dominate over the linear growth factor leading to a net increase over the 
case with the same Qm- However the enhancement due to the integral over x has to contend 
with the factor of outside the integral which is considerably larger. The net result is 
shown in figure |1] for source galaxies at Zg = 1 (upper panel) and Zg = 3 (lower panel). The 
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y-axis shows times the x-integrand of equation (||). Thus Cpp or is proportional to 
the area under the curves shown, aside from numerical factors which depend on the shape 
of the spectrum. The figure shows that the peak contribution for a reasonable choice of Zg 
comes from the range z ~ 0.3 — 0.6. 

The scaling of p{6) (the same as that of the square root of Cpp) with VL^ is dependent 
both on Zg and the shape of the power spectrum. For Vtf^ = the dependence is well fit by 
the power law, 

Zs = l : p{e) cc asnl''. (10) 

For Zs — I and — 2 n ^ — 1, the above equation is quite accurate. For Zg = 3 the scaling 
depends on the spectrum, and can be approximated as, 

Zs = 3 : p{9) oc agl^J^^^ for n = -I ■ p{e) oc ag^^^^^ for n = -2 . (11) 

The result for n = —2 is valid provided a low-Zc cutoff is imposed to keep the integral finite. 
As expected, at higher source redshifts the enhancement due to the distance and growth 
factors is larger and leads to a weaker flrn dependence. The same occurs as n decreases, 
as then there is more power on large scales, and therefore more weight to the high- 2; part 
of the integral. The results of equations ( p!0|) and (0) are valid for 0.2 ^ Qm ~ 1, and are 
intended to cover the range of source redshift and spectra spanned by realistic models. 

A simple way to understand the dependence on Qm follows from the fact that most of 
the lensing contribution comes from about half the distance to the source galaxies. Thus the 
flm dependence of equations (p!0| ) and (|Tl]) can be compared to that of Qm at the redshift 
Zi/2, corresponding to Xsf^- ^m(-2i/2) is well approximated by fl'^ for Zs = I and by fl^ 
for = 3 if n = —2. For a shallower slope there is more power on small scales, which are 
closer for a given angle. Hence the results are closer to those of equations (|10|) and ([TT|) for 
n = — 1 if one uses Xs/S (it makes very little difference if one uses r{x) instead of x for this 
purpose). Thus as a first approximation the flm dependence can be obtained by replacing 
the Qrn term in Cpp for the flat case, by the value of at 1/3 to 1/2 the distance to the 
source galaxies. 

The above scalings are derived using linear evolution of the power spectrum. On smaller 
angular scales nonlinear evolution could play an important role — in the next section 
we compute the change in the above scalings due to nonlinear effects. It is complicated 
to estimate the angular scale where nonlinear effects are important because each angle 
involves the projection of a range of length scales at different redshifts. The answer will 
therefore depend on the slope of the power spectrum. A rough estimate of the length 
scales that make the dominant contribution to p{6) is given by figure 0. This figure shows 
the transverse distance for a given 6 at the redshifts estimated from figures |^a and 0b as 
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providing the dominant contribution to p for source galaxies at ~ 1 — 2. For ^ < 10' this 
length scale is less than 5h~^Mpc, which is the correlation length of galaxies and therefore 
a good demarcation of the nonlinear regime. We may thus expect that the linear results for 
p will be modified for 6 < 10'. 

An analytic estimate of the change in the above scalings due to strongly nonlinear 
effects can be made using the stable clustering solution. If one assumes that the lensing is 
dominated by clustering on small scales (valid in the limit of small angles and low source 
rcdshift), which stabilized at high redshift. then the dimcnsionlcss power is a function of the 
physical wavcnumbcr times a^. This property extends to the open models as well, because 
at high enough redshift, flm{z) ~ 1. It can then be shown that, if the linear spectrum is 
normalized to the same erg today, then the enhancement factor of the nonlinear spectrum 
relative to the linear one is proportional to F{x)~^- This follows from a reversed version of 
the argument of Peacock & Dodds (1996), which was made for spectra normahzed to the 
same initial value. A particularly simple example is provided by the n — —2 spectrum, for 
which the nonlinear slope and growth in the Einstein-de Sitter case are the same as in the 
linear regime. Hence the nonlinear spectrum is given at all /c by a constant enhancement 
factor relative to the linear spectrum. 

For 0.3 ^ flm ~ 1 one can approximate F{x) — for — 1- Therefore the 
Qm dependence in p'^{6) due to strongly nonlinear evolution of the power spectrum is 
-^(x) ~ ^m^'^- ^^^^ regime then, p"^ oc fi^^ times the contribution of the distance 
factors, which further lowers the power of Qm- This estimate is meant to provide an 
upper bound on the enhancement for open models due to the nonlinear contribution. For 
reahstic redshifts the contribution from such small scales does not dominate, and the exact, 
numerical estimates of the nonlinear spectrum must be integrated over k and x- Still, it is 
clear that the D,rn dependence of weak lensing amplitude will be significantly weaker than 
the linear relation predicted from the scalings in previous subsection. 



3. Predictions for rms ellipticity for CDM-like spectra 

The detailed predictions for the rms ellipticity depend on the shape and amplitude of 
the power spectrum, the distribution of source galaxies, and on the cosmological parameters 

and r^A- In the following subsection we shall consider two alternative normalizations of 
realistic CDM-like power spectra. A useful, model independent description of the scalings 
can also be obtained by fitting the dependence on the different parameters to power laws. 
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Table 1: Scaling of the polarization p{6) with 6, Zs and with cosmological model parameters 
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This provides a reasonably accurate approximation to the result, and shows the qualitative 
behavior more clearly. 

Table 1 provides such power law fits to the scalings of p with 6, Zg, cg, and Qm for 
Einstein-de Sitter, open and A— dominated cosmologies. The three specific models we 
have chosen are: Einstein-de Sitter [Q^ = 1? = 0), Open [Q^ = 0.3, Q\ = 0), and 
A— dominated {Qm = 0.3, i^A = 0.7). The latter two models are representative of the class 
of open and fiat A— dominated cosmologies as the results do not change significantly for Qm 
in the range 0.2 < < 0.5. The matter power spectrum used in the calculations is the 
nonlinear F = 0.25 CDM spectrum, where F is the shape parameter defined in Bardeen et 
al. (1986) and roughly corresponds to Qmh in CDM models. The scalings given in table 
1 are not strongly sensitive to the shape of the power spectrum as shown in figure |^ and 
discussed in the next subsection. 

On scales 6 < 10' nonlinear evolution makes a significant difference to all the scalings. 
Figure || shows the ratio of p{9) computed using the nonlinear/linear power spectrum for 
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the above three models. The ratio is significantly larger than unity for 9 < 10', and reaches 
a factor of two for 9 ~ 1'. It is the largest for the open model, where the nonlinear effects 
are more important because of a slower linear growth of perturbations. 

Since the scalings differ in the linear and nonlinear regimes, table 1 gives power law 
fits for three different angles, 9 — 2', 5', 15', to represent the nonhnear, weakly nonlinear 
and linear regimes, respectively. These give p{9) to good accuracy over most of the range: 
V < 9 < 30'; 0.5 < < 3; 0.5 < cig < 2; 0.2 < Q^n < 1- The results can be approximated by 
the following equations for the three cosmological models given above. The angle 9 below 
is in arcminutes, and where the scalings differ on 2' and 15', the latter scaling is given in 
brackets in italics. Since the scaling with Zg is nearly the same for the two values of 6*, we 
have used an intermediate value in the following equations. 

P[2'ii5')] = 0.08 r°-=^^ ag^-^'^') : nr^ = i.o 

P[2'{15')] = 0.04 r"-^^ 6 ag'-^'^'^ (l]„/0.3)0-*^6(''-^^) : 17^ = 0.3 
P[2'{15')] = 0.04 r°-^2 ^0-^6 ^i.29(i) ^^^^Q_3^o.60(o.es) . Oa = 0.7. (12) 

These power law fits can be used to obtain predictions for the rms ellipticity for any 
desired choice of cosmological parameters. For 9 — 15', the results are consistent with the 
linear theory estimates for power law spectra (with n between —1 and —2) obtained in the 
previous section. The p{9) curve is enhanced at small angles due to nonlinear evolution. 
This enhancement causes the dependence on 9 to be much closer to a power law over the 
range V < 9 < 30' than the linear prediction, which has much greater curvature due to the 
curvature of the power spectrum. This power law is steeper for the open and A models due 
to stronger nonlinear evolution. 

On 9 ~ 2', the dependence on as is significantly stronger than linear, especially for 
the open model. Nonlinear evolution also weakens the dependence on Qm as anticipated 
in the previous section. Combining the scaling with <Ts and Qm (relevant for comparison 
with other tracers), p{9) for ~ 2' scales as (Xg'^^ ^m^^j whereas on large scales of > 10', it 
scales as erg Thus on small scales p{9) measures crgfl^ (the power of flm is even lower 
for ^ < 2' or Qm < 0.5), whereas on large scales it measures ag 

Using the above differences a comparison of p{9) at small and large angular scales 
can constrain ag and Qm separately. Thus nonlinear evolution, at the expense of a more 
complex dependence on various parameters, can help break the degeneracy between these 
parameters. To distinguish an open from a flat A— model with the same value of one 
needs to use the difference in the scaling with Zs, as the other scalings are very similar. The 
A— model predicts a faster increase with Zg due to a faster increase of the distance factors 
with redshift than the open case. Once it becomes feasible to get estimates of the redshifts 
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of source galaxies, this difference can become a useful test of A as it is nearly independent 
of the shape and normalization of the power spectrum. 



3.1. Results for cluster-abundance and COBE normalized spectra 

In choosing the most realistic power spectrum to compute the polarization signal, the 
approach closest to other tracers of the mass density field would be to use the abundance 
of galaxy clusters or the peculiar velocity field of galaxies to estimate the amplitude erg, 
and to then use the galaxy power spectrum to constrain the shape (Peacock & Dodds 
1996; Peacock 1996). The merit of this approach is that it agrees with the observational 
constraints on the scales which are most relevant for weak lensing observations. An 
alternative is to work with physical models for the dark matter and cosmology which give 
the shape of the power spectrum, and use COBE-normalization to fix the spectrum on 
very large scales and thus its overall amplitude. The merit of this approach is that for a 
given model both the spectrum of perturbations and their evolution is fully determined. 
However, the final spectrum may or may not agree with the present day observations (which 
interpretation is still somewhat uncertain). We shall present results using both approaches. 

The first of the two reliable tracers of mass fiuctuations on large scales is the abundance 
of large galaxy clusters with a typical mass of lO^^M©, which corresponds to a linear scale 
of about 8/i~^Mpc. These objects are very rare and according to the Press- Schechter 
formalism (Press & Schechter 1974) their number density depends exponentially on the 
amplitude of mass fluctuation on the corresponding linear scale. This allows an accurate 
estimate of the amplitude, which typically gives as — 0.5 — 0.6 (White, Efstathiou & Prenk 
1993, Viana & Liddle 1995, Eke, Cole & Frenk 1996, Pen 1996) in an = 1 universe, 
almost independent of the shape of the power spectrum. At present this gives the best 
constraint on the linear density power spectrum, with an error of about 10-20%. The 
dependence on is approximately cTg oc and differs only slightly for open and 

spatially flat models. The scaling with flm and fl\ is still somewhat uncertain as it requires 
calibration with large N-body simulations. 

Peculiar velocities are also tracers of underlying mass fluctuations. In the 
linear regime there is a simple relation between velocity and density fields (Peebles 
1980):5(r) = — (ifo/)~^V • v{r), where / ~ fl^. The Vim dependence on mass fluctuation 
amplitude from velocity data is as oc which is almost the same as for the cluster 

abundance in a flat universe and only slightly less so in an open universe. Comparison of 
the two constraints therefore cannot signiflcantly constrain Q^. Given a set of velocity 
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measurements one can reconstruct the particular combination fS and estimate its power 
spectrum. Recent comparisons between Mark III TuUy-Fisher catalogs and IRAS survey 
tend to favor pj = n^-^bi ~ 0.5 ( |Willick et al. 1996| , Pavis, Nusser fc Willick m^ ), where 
bi is the linear bias parameter for IRAS galaxies. This combined with ag.i ~ 0.7 leads to 
agfl^'^ = 0.35, smaller than the value derived from the cluster abundance data. This would 
therefore reduce the weak lensing predictions presented in this paper. However, analysis of 
Mark III catalog with POTENT gives a higher value of as^°-^ (Kolatt & Dekel 1996), so it 
seems prudent to adopt the cluster abundance normalization until the discrepancy in the 
velocity data is resolved. 

Based on the above discussion, a useful functional form for the density power spectrum 
is given by a CDM type transfer function with two free parameters, the amplitude 
as ~ 0.6^"°'^ and the shape parameter F ^ 0.25 (Peacock & Dodds 1996). This is 
the model that will be adopted in computing predictions of ellipticity polarizations for 
cluster-abundance normalization. We shall use the fitting formulae in equations (47-49) of 
Viana fc Liddle (1995)| to compute (Xg as a function of Qm ^ their normalization is very close 



to the results of White, Efstathiou fc Frenk (1993) and Pen (1996), though a bit higher 



than that of |Eke, Cole fc Frenk (1996^ 



Before proceeding with the results for p{6) we consider the contribution to the integral 
from different ranges in k. Figure § compares the contribution per logarithmic interval in 
wavenumber to and cr^(/) at ^ = 15', Zg = 1 with the logarithmic contribution to cr|, which 
is also an integral of the power spectrum over k with a real space top-hat filter. All the 
distributions are normalized so that they peak at unity. One can see that the distribution 
of is broader and that of cr^(Z) narrower than the corresponding ct| distribution, but in 
general they sample very much the same scales. Thus the optimal statistic to use in order 
to compare the normalization as with other tracers is cr^(/), p'^{0) at 1//, ~ 15'. The 
comparison can then be made nearly independent of the shape of the power spectrum. 

We need to choose a redshift distribution for the source galaxies in order to compute 
p{6). At present this is spectroscopically known only for galaxies with magnitude in I 
below 23 (Lilly et al. 1996; Cowie et al. 1996), which have median redshift z ^ 0.5, while 
typical lensing observations reach several magnitudes deeper (e.g. Mould et al. 1994). 
HST observations of Abell 2218 coupled with cluster mass reconstruction indicate that 
the median redshift increases from z ~ 0.5 at R=22 to 2; ~ 1 at R=25 (Kneib et al. 
1995). There could however be a population even at redshift beyond 1, as indicated by the 
detection of the shear around a z = 0.8 cluster (Luppino & Kaiser 1996). Most of the results 
we present are for Zg = I, but they can be adapted to any desired Zs using the scalings of 
table 1. For simplicity we have assumed that all the galaxies are at the same distance. This 
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is the best possible case in the sense that it makes the distributions in k space more narrow 
than for more reahstic cases. Fortunately the difference between our results and those using 
a different, more realistic distribution of source galaxies with the same median redshift is 
usually very small (Kaiser 1992, Villumsen 1996). 

Figure ^ shows the dependence of p{d) on 6 for the three cosmological models, all 
with shape parameter F = 0.25, except for the Qm = 1 COBE-normalized model. The 
thick solid curve is for Q^n = 1? the dashed curve for Q^n = 0.3 and the dotted curve for 
Qm = 0.3, Qa = 0.7. The normalization chosen is (a) as = 1, (b) Cluster-abundance, and 
(c) COBE. For the open and flat COBE models, we have used Qm = 0.4 so that the value 
of h and erg are reasonable. The three COBE models are: flat, = 1, tilted with n = 0.8, 
h = 0.5, as = 0.72; open with fi„ = 0.4, h = 0.65, as = 0.64; fiat, A-model with = 0.4, 
h = 0.65, as = 1.07 (Bunn & White 1996). Since the A— model chosen has a high as, it 
turns out to have nearly the same p{9) prediction as the tilted Qm = 1 niodel. 

The thin solid curve is computed using the linear spectrum for the Qm = 1 niodel. As 
expected from the results of figure |^, it is only for 9 > 10' that the linear and nonlinear 
curves start to coincide. For cluster-abundance normalization, the amplitude of ellipticity 
correlations very weakly depends on the mean density for a reasonable range of Qm- The 
enhancement for low Qm models from the growth and distance factors discussed in §2 nearly 
cancels out the fl'^ dependence one might have expected from the term outside the integral 
in equation |^. Therefore one has to make subtler comparisons of the signal at small and 
large 9 to break the degeneracy between as and flm- This was discussed above following 
equation ([T^). It indicates that comparing weak lensing measurements over different scales 
can more powerfully constrain flm than comparing the weak lensing amplitude with other 
tracers. For COBE normalization on the other hand, a detection of the correlated ellipticity 
with reasonable error bars can immediately constrain the cosmological model in question, 
but of course there might still remain several different models which equally well fit the 
data. 

The effect of choice of statistic, shape of the power spectrum, and the redshift of source 
galaxies is shown in figures ^ and ^ Figure ^ shows a{l/9) for the three cosmological 
models with cluster-abundance and COBE normalization as in figure ^. The curves for a 
flatten at small 9 more than the p curves. This is to be expected because p integrates over 
all scales larger than r9 (where r is the typical distance) and so can only increase toward 
small scales. Because there is little power on very small scales it eventually saturates (this 
typically happens on subarcminute scales, Keeton, Kochanek & Seljak 1996). The Fourier 
space quantity cr^(/) on the other hand receives contributions only from scales smaller than 
r/l and therefore decreases on very small scales for realistic power spectra. The shape of 
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curve is very similar to that of p^, and for 9 2', Cpp is well approximated by 0.7 p^. 
Figure |^ compares p{6) for the F = 0.5 CDM spectrum to the F = 0.25 spectrum used in 
the rest of the figures. The upper and lower set of curves are for Zg = 2 and 1 respectively. 
Since the F = 0.5 spectrum has more power at small scales, it leads to a larger p at small 6. 
The curves cross around 6 = 10', but the difference is not large for 6 > 2'. Figure || shows 
the variation of p{6) with for 6 = 15' and 6 = 2'. The models shown are the same as in 
the cluster-abundance normalized panel of Figure ^. As shown in Table 1, p increases with 
Zs fastest for the A— models. 



4. Nonlinear effects on the distribution of ellipticity 

So far we have only discussed the importance of nonlinear effects on the second 
order statistic, which as we have seen increases the signal above the values predicted by 
linear calculations. Nonlinear evolution however also changes the distribution function 
of correlated ellipticity or magnification. In particular, the tails of distribution grow 
and become asymmetric, reflecting the nonlinear growth of perturbations into small and 
overdense objects. This leads to a distribution of magnification which is skewed toward 
positive values. Such skewness is of interest for determining cosmological parameters by 
itself (Bernardeau et al. 1996), but also affects the inference of cosmological parameters 
from the measurements of 2-point statistics. If the distribution of polarization at a given 
angular scale 6 can be assumed to be normal distributed then the usual results from 
gaussian statistics apply. For example a single measurement of polarization amplitude 
within a circular aperture has a Rayleigh probability distribution x^(2), which has 
exponentially suppressed tail and from which the limits on cosmological parameters can be 
derived analytically. Several independent measurements can be easily combined and only 
a few observations suffice to obtain strong limits on cosmological parameters. If, on the 
other hand, the distribution is strongly non-gaussian so that the tails of distribution are 
important then one cannot use these simple arguments to derive the limits and a larger 
set of measurements (a "fair sample") is needed to sufficiently sample the tails so that the 
variance converges to a true value. At the same time, higher order statistics become easier 
to measure and can provide further tests of cosmological parameters. 

For measurements of density perturbations a simple criterion of where the nonlinear 
effects become important is cij? ~ 1, so for scales above R ~ 10h~^ Mpc the universe is in 
the linear regime and gaussian statistic can be used. For smaller scales nonlinear effects 
become important which leads to deviations from a normal distribution even if the initial 
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field was gaussian. For projected quantities such as shear and magnification the answer 
is not so simple. The projection integrates over all the scales, so there will always be 
contribution from very small and thus very nonlinear scales. On the other hand, projection 
will reduce the deviations from a gaussian distribution because of the central limit theorem: 
if we sum over sufficiently large number of deflectors the flnal distribution will be a gaussian 
regardless of the probability distribution of individual deflectors. We thus expect that the 
nonlinear effects will be more important for small angles and for bright (nearby) background 
galaxies. 

To obtain a more quantitative answer we calculated the skewness of local convergence 

5 = using second order perturbation theory. Skewness for local convergence 
does not vanish as in the case of ellipticity (Bernardeau et al. 1996). The second moments 
of 2R, and p are the same and so one may expect the higher order moments of convergence 
to be indicative of the nonlinear effects on the ellipticity distribution as well. Note that the 
quantity S differs from the quantity S3 = (k^) / (S^)^ used by Bernardeau et al. 1996. The 
latter is deflned so that it is independent of the amplitude of power spectrum in the second 
order perturbation theory, while S* is a more useful variable to quantify the deviations from 
the gaussianity. As discussed in the appendix, second order perturbation theory should 
be approximately valid in the regime of interest here, between 1' to 1°. Second order 
perturbation theory accurately estimates the skewness over a much wider range of scales 
than it does the variance. Only on the smallest scales it underestimates the skewness, but 
typically not more than a factor of two, as shown by N-body simulations (Colombi, Bouchet 

6 Hernquist 1996). Using the expressions derived in appendix (see also Bernardeau et al. 
1996) we can compute the value of S for any of the models discussed in this paper. Figure 
P presents the distribution of as a function of 6 for 3 different redshifts in 3 different 
models: flat model with ag = 0.6, curvature and cosmological constant models, both with 
CTg = 1 and fl = 0.3 (this is very close to the cluster-abundance normalization). In all the 
models we use F = 0.25. 

For simplicity we will choose S = 0.2 as the value where the distribution becomes 
non-gaussian. For background galaxies at Zg = 0.5 we flnd that the distribution is 
non-gaussian even at 4° in the flat model, where the polarization signal is only 0.25%. At 
Zs = 1 and Zs = 2 the corresponding values are 1° and 0.5°. For curvature and cosmological 
constant dominated models the non-gaussian effects are even more important because of 
higher normalization at Sh^^Mpc, although the effect is offset by the longer radial distance 
to a given redshift. The conclusion is that in most of the observationally interesting 
regime the nongaussian effects will be quite important. While this complicates the issue 
of how large a fair sample should be, it also implies that skewness will be fairly easy to 
measure with a reasonably small sample. Because S3 depends sensitively on Qm but not 
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on the amplitude of the power spectrum (Bernardeau et al. 1996) this test may become a 
promising way to determine Qm- 

The analytical approach presented above provides only a limited view of the properties 
of ellipticity distribution. For example, such an analysis cannot answer whether skewness is 
produced by a few outlying points generated by rare clusters, in which case a rather large 
sample would be needed to sample these properly. For a quantitative estimate of what 
constitutes a fair sample for measuring the variance requires therefore knowledge of at least 
the fourth moment. Such estimates of a fair sample for measuring the second and third 
moment can be made more directly by examining the distribution of correlated ellipticities 
in N-body simulations (Wambsganss, Cen & Ostriker 1996). 



5. Discussion and conclusions 



Weak gravitational lensing has several advantages over other probes of large scale 
structure. It is directly sensitive to the underlying mass distribution and is therefore 
independent of biasing, which has plagued many of the cosmological tests based on the 
distribution of galaxies. While peculiar velocities in principle similarly probe the dark 
matter distribution directly, in practice they suffer because of large observational errors and 
various biases associated with these. This is the main reason why higher order statistics of 
velocity fields have not provided strong constraints on cosmological models. 

Gravitational lensing could provide a much cleaner way of performing power spectrum 
and higher order statistical tests, provided that one can overcome observational difficulties 
and reach the 1% level of polarization, quite feasible with the new generation of composite 
CCD cameras. For a filled survey with an area 6"^ and measured galaxies the noise 
variance in p'^{6) is given by ~ (0.4)^6'o/A^^, where 0.4 is the typical intrinsic shear for a 
single galaxy. Although one has to average over sufficient number of galaxies to reach this 
level of signal, this is not a severe limitation for sufficiently deep exposures with several 
hundred thousand galaxies per square degree, so on degree scales the signal to noise is of 
the order of 100 (Kaiser 1996; fig. |l^) and in fact sampling variance (finite number of 
independent areas of size 6"^) is more important than noise on large angles. On smaller 
scales the noise increases (proportional to 6~^), but the signal also increases nearly as 
rapidly, so that the overall signal to noise remains approximately fiat down to 2' scales in 



open and cosmological constant models and slowly decreases in fiat models (fig. |T0]). Based 
on this it appears there should be sufficient power present in a survey of this size to measure 
the signal even on arcminute scales with a high statistical accuracy, so a comparison 
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between small and large scales to constrain would indeed be feasible, provided that the 
systematic effects can be kept under control. 

The other limitation at present is that the redshift distribution of background galaxies 
is unknown beyond / = 23 (Lilly et al. 1996; Cowie et al. 1996). New observations with 
Keck will be able to extend this spectroscopically up to / = 25 (e.g. Lowenthal et al. 1996) 
and even beyond, if spectra of magnified arcs are measured. In addition, weak lensing on 
high redshift clusters can be used to constrain the redshift distribution of faint galaxies 
without measuring their redshifts, although at present there appears to be some controversy 
over whether the faintest galaxies are predominantly around or below z = 1 (Kneib et al. 
1995) or above that (Luppino & Kaiser 1996). Finally, the best hope to constrain their 
redshift distribution may lie in the photometric redshifts technique (ConoUy et al. 1995; 
Sawicki, Lin & Yee 1996). This would allow one to perform such precision tests as the 
growth of power spectrum with redshift, needed to estimate the value of ^a. For example, 
based on Hubble Deep Field observations, Sawicki et al. 1996 find that the redshift 
distribution of / < 27 galaxies is bimodal, with a large contribution of (predominantly 
blue) galaxies between redshifts 2-3, while at / < 25.5 most of the galaxies are below 
z = 1. In both cases the weak lensing signal to noise with these limiting magnitudes would 



be comparable to the one in fig. |10|, so in principle by comparing the signals in the two 
magnitude bins it would be straightforward to deduce the value of ^a. The caveat is that 
the faintest galaxies in HDF appear to be very irregular and small, so it is not clear whether 
weak lensing analysis from the ground could be performed on these objects. While more 
investigation is needed along these lines, there appears to be a real promise of measuring 
cosmological parameters with future weak lensing observations. 

We have analyzed theoretical predictions for weak lensing ellipticity correlations over a 
broad range of angles and for a variety of cosmological models. Our principal conclusions 
are the following: 

a) the amplitude of the rms ellipticity varies between 4-6% on arcminute scales to 
0.5-1% on degree scales for background galaxies at Zg = I and grows as In this 
range the amplitude only weakly depends on the value of cosmological parameters if cluster 
abundance (or velocity flows) normalization is adopted. A single measurement of weak 
lensing amplitude on ^ ^ 15' will provide a determination of erg ^^m~^'^ which can be 
compared with the cluster-abundance or velocity flow predictions. Because these tracers 
all scale similarly with this will not provide an independent determination of flrn, for 
which a more complicated analysis will be needed. 

b) nonlinear effects are important for second order statistics on scales below 10' and 
enhance the signal expected from linear theory. This enhancement can be up to a factor of 
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two on arcminute scales. From this it follows that the signal to noise does not decrease as 
rapidly on small scales as expected from linear theory and is in fact approximately flat in 
certain models from arcminute to degree scales (see figure |T^). Therefore measurements over 
this whole angular range should be pursued, especially since they provide complementary 
information. The interpretation of cosmological models is more cleanly addressed on 
large angular scales, where nonlinear effects are negligible and the distribution function 
is primordial. Small angular scales can provide insight on the nonlinear evolution of 
perturbations. If our current understanding of the latter is correct then combining 
measurements on large and small angular scales will help break the degeneracy between the 
cosmological parameters, because of different scalings in the linear and nonlinear regimes. 
The ellipticity on small angles measures agfl^, on large angles it measures asfl'^, while 
the variation with Zs can constrain f^A. Thus, to obtain model independent measurements 
of the cosmological parameters it is necessary to measure the signal over a range of angles 
and source redshifts. 

c) nonlinear effects on the distribution of ellipticity are significant on scales below 1° 
and imply that gaussian statistics cannot be applied to the data. This means that a larger 
sample will be required to measure second order statistics than expected from a gaussian 
distribution. Conversely, skewness and other higher moments will be easy to measure and 
may provide a strong test of cosmological parameters independent of the tests above. 

In conclusion, weak lensing has the promise to become the next testing ground for 
cosmological theories. The expected signal to noise is above unity over a large range of 
angles and limiting magnitudes. This should allow one to test sensitively the spectrum 
of density perturbations and its evolution in time, in both the linear and the nonlinear 
regimes. 

B.J. is very grateful to Jens Villumsen for initiating his interest in weak lensing, and 
to him and Richhild Moessner for useful discussions. We thank Bernard Fort, Nick Kaiser, 
Chris Kochanek, Yannick Mellier, Stella Seitz, Simon White, and especially Peter Schneider 
for stimulating discussions. We thank the referee for suggestions which helped improve the 
paper and led to the addition of figure 8. 



A. Appendix 

In this Appendix we present a general description of lensing in a weakly perturbed 
universe. The derivation presented here follows Seljak (1995, 1996) and differs from recent 
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derivations in Bernardeau et al. (1996) and Kaiser (1996) in that it employs a global 
description of photon trajectories. The geometrical interpretation requires only knowledge 
of spherical trigonometry and it is thus particularly simple to derive the general weak lensing 
expressions. Moreover, it can easily be generalized beyond weak lensing approximation or 
small deflection angles. 

The framework is a perturbed Robertson- Walker with small-amplitude metric 
fluctuations. We will only consider scalar perturbations here, in which case the metric in 
the longitudinal gauge can be written using conformal time r and comoving coordinates x'^ 
as 



Here a(r) is the scale factor expressed in terms of conformal time. We will adopt units such 
that c = 1. We set the two scalar potentials to be equal, which is a good approximation in 
the matter dominated epoch. The space part of the background metric can be written as 



where K is the curvature term which can be expressed using the present density parameter 
Vt (only the present day density parameter will be used in this paper, hence we may drop 
the subscript 0) and the present Hubble parameter Hq as K = {VL — 1)Hq = —QkHq. 
The relation between the radial distance x ^-nd the redshift z can be obtained from the 
Friedmann equation, da/dx = HQiVtrnO^ + Vtf^a^ + fi^^a^)^/^, where a = (1 + z)~^ is the 
expansion factor and the densities of matter, cosmological constant and curvature are 
expressed in terms of the critical density. The density parameter Vt can have contributions 
from mass density VLm or vacuum energy density ^a, ^ = + ^a- The advantage of using 
the conformal time r is that the metric becomes conformally Euclidean {K = 0), 3-sphere 
{K > 0) or 3-hyperboloid {K < 0) and leads to a simple geometrical description of light 
propagation. 

In an unperturbed universe a photon emitted from a source toward the observer will 
travel along a null geodesic in the radial direction with a radial position given by x = — t- 
Adding a perturbation changes the photon trajectory. The change in photon direction is 
governed by the space part of the geodesic equation, which applied to the metric (equation 
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-fijdx'dx^ = dx^ + r\d9^ + sin^ Odcp'^), 



r(x) = sin/^x 
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where the symbol V_l0 denotes the transverse derivative of potential. Null geodesies obey 
ds^ = 0, from which follows the relation (ix = (1 — 2(f))dl. Here can be interpreted as 
the Newtonian potential, since on scales smaller than the horizon it obeys the cosmological 
Poisson equation. It can be viewed as providing a force deflecting the photons and affecting 
their travel time while they propagate through the unperturbed space-time. Equation ( |A^ ) 
is a generalization of Einstein's deflection angle formula and includes the well-known factor 
of 2 difference compared to Newtonian gravity. Even when metric perturbations are present, 
one can continue to parameterize the geodesic with the unperturbed comoving radial 
distance x- The deflection angle at a given position x can be calculated using a locally fiat 
coordinate system, which allows a plane wave expansion for the potential 0, provided that 
the longest correlation length is small compared to the curvature length (this condition is 
well satisfied for the power spectrum in our universe). The effect of the deflection angle on 
the photon transverse position must however include the curvature effects. In practice this 
means that one only needs to know how to solve triangles using the spherical, Euclidean or 



hyperboloid trigonometry (figure |TT]). Because the only observable photon direction is that 
at the observer's position we will propagate photons relative to their final direction (i.e. 
backwards in time). We will also adopt a small deflection angle approximation, because 
one does not expect large deflection angles due to the lensing, but the expressions can be 
generalized to remove this restriction. In this approximation the transverse derivatives in 
equation ( [A3| ) can be approximated with the transverse derivatives with respect to the 
observed direction of the photon or with respect to any other fiducial direction that has a 
small angular separation with the photon (for example the unperturbed direction to the 
source). In this plane approximation the observed photon direction n can be described with 
a two-dimensional angle 6 with respect to the fiducial direction, n = {6i, 62, 1 — |^P/2 ~ 1). 

Suppose a photon is observed at an angle 9 relative to some fiducial position. As it 
propagates through the universe the photon is additionally deflected according to equation 
(see figure |1T]). This leads to the transverse photon excursion x±{x) relative to the 



unperturbed line. From Euclidean, spherical or hyperboloid trigonometry one finds that an 
individual deflection by 6d at x' leads to an excursion at x given by 6x± = Sdr{x — x')- 
The total excursion is given by an integral over individual deflections, 

x^ix) = -2 r V^(j){x')r{x - X)dx + Hx)- (A4) 

JO 

Because of gravitational lensing the "true" surface brightness (i.e. the surface 
brightness one would see in the absence of any lensing) at position xs is mapped into the 
observed one, Iohs0) = hrue{d + S9), where 59 is the angular deflection of a photon caused 
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by intervening mass, which follows from equation (|A4|) 

Se = ^^ rixs~x)^±<Pix)dx- (A5) 
r{Xs) Jo 

While the deflection 69 is not directly observable, its gradient is, through the stretching and 
magnification of distant galaxies. This is described by the two-dimensional shear tensor 
defined in equation (|l]). To compute it we expand the potential in equation (|A5|) across 
two neighboring rays separated at the observer by d6 using the unperturbed separation 
dx = r{x)d6. This approximation assumes that the components of shear tensor are 
small (Kaiser 1992), or, similarly, that the relative deflection between the neighboring 
rays is small compared to the unperturbed separation. This is the so-called weak lensing 
approximation. 

The shear tensor can be diagonalized and decomposed into its trace 2k = $ii -|- $22 
and ellipticity (or polarization) p = pi + ip2 given by 

p^ = - $22 ; P2 = 2$i2 . (A6) 

While trace 2k magnifies (or demagnifies) the images of galaxies, ellipticity p stretches the 
images and can be observed by averaging over a sufficient number of galaxies so that the 
noise caused by the intrinsic ellipticity of galaxies becomes smaller than the signal. The 
simplest average is the one within a circular window of radius 6, which can be written in 
Fourier space as 

PiO) = d'9'J^ g{x)dx J dV-^e'''--''klf3m. (A7) 

We decomposed the wavevector into the radial component and the 2-dimensional 
transverse component k± with the amplitude k_i_. We introduced the variables 
f3 = exp{2i(f)k), where (pk is the azimuthal angle of k±. The integral over d'^6' can be readily 
performed to give 

rX.0 r 

p(9) = 2J^ g{x)dx J d^ke'''^''klf3(l){k)W2{k^r9), (A8) 

where W2{x) = 2Ji{x)/x. We may now employ the small angle approximation (Limber 
1956), which is valid in the limit where the radial window function is broad compared 
to the typical wavelength that contributes to the integral. In this limit only the modes 
perpendicular to the radial direction will contribute to the integral, all the others being 
suppressed because of cancellation of positive and negative fluctuations along the line 
of sight. This allows one to put k± ^ k. With this approximation one obtains for rms 
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ellipticity amplitude p = pp 



{p\d)) = 4 / / d%d% 
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(A9) 

Ensemble average can be expressed in terms of power spectrum of potential P,p{k) or density 
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where power spectrum P5{k, a) and expansion factor a depend on time. Integrating over the 
Dirac do function and over the radial component of wavevector we finally obtain equation 
(0), which gives the rms amplitude of polarization. The same expression can also be used 
for rms amplitude of magnification 2R. The rms amplitude for each of the two ellipticity 
components is simply 1 / \/2 of the total. 

Because skewness of ellipticity vanishes for symmetry reasons we will use skewness 
of local convergence as an estimate of the nongaussian nature of ellipticity distribution. 
Calculating it proceeds along the same lines as for the variance (Bernardeau 1995). It is 
defined as an ensemble average of the third moment of the mean magnification, 



{>e{e)) 



xo 



dx / d^ke'''^''g{x)W2{k^r9)k^4>{k) 



(All) 



In linear theory the 3-point correlation function vanishes, so one needs to go beyond the 
linear approximation to obtain a nonvanishing value. In second order perturbation theory 
the density field is given by 



^(2)(fc) 
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d%d%S^^\ki)S^^\k2)SDih + k2- k)F{h, h 
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Note that the time dependence is all in the growth factors of 5{k)^ which is calculated using 
linear theory (neglecting the extremely weak dependence of F{ki,k2) on Q, Bouchet et al. 
1992). Ensemble averaging of skewness gives 6 identical terms, which after integrating over 
the radial wavevectors leads to the following expression, 

{R^e)) = Sln'e-'H', r g^D'n'a-'r-^dx 

Jo 

' d%d%W2{li)W2{l2)W2{\k + l2\)P6{klre)Ps{l2/re)F{lj2). (A13) 
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The integrals can be further simphfied using similar manipulations as in Bernardeau (1995), 
which leads to the expression (see also Bernardeau et al. 1996) 

{k%9)) = ?,22Tx^d-^Hl j^\^D^(^\-^dx j P{l/re)Wi{l)ldl 

(AM) 
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^ / P{l/re)W^{l)ldl + i / P{l/re)W2{l)W!i{l)l'^dl 
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This expression was used in §4 to evaluate the importance of nonlinear evolution on the 
distribution function of polarization. 

Several approximations have been employed in obtaining the result above, which limit 
its validity. Second order perturbation theory results are only valid in the domain where 
density perturbations are not much larger than unity, so on very small scales a more 
involved calculation would be needed. However, these corrections are typically not more 
than a factor of two. Because our purpose is to estimate at which scale the nongaussian 
effects become important second order perturbation theory suffices. Another approximation 
is the use of small angle approximation, which leads to corrections on angular scales above 
1° (Bernardeau 1995). In the regime of main interest for us, which is between 1' — 1° this 
is not an important correction. Finally, expression (|A14[ ) neglects the deviation of the 



photon trajectory from the unperturbed path, which also contributes at the second order. 
This correction is less important than the main term used in equation (|A14D and can be 
neglected (Bernardeau et al. 1996). 
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Fig. 1. — The dependence of the integrand for Cpp on redshift. The source galaxies are 
assumed to be at 2;^ = 1 (upper panel), and Zg = 3 (lower panel), and the matter power 
spectrum is a power law with slope n = —2. The solid curve is for Qm = Ij the dashed 
one for Qm = 0.3, Q\ = 0, and the dotted one for Qm = 0.3, = 0.7. The figure shows 
that compared to the Einstcin-dc Sitter case, the amplitude decreases for the A model and 
is lowest for the open model. The peak of the integrand for Zs = I lies at about z = 0.3 
for the — 1 case, and shifts to higher z for the open and A models. The light, dashed 
curve is for an n = — 1 spectrum, with = 0.3. It peaks at lower z compared to the heavy, 
dashed curve as the spectrum has more small scale power. 
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Fig. 2. — The transverse distance at z = 0.3 (lower set of curves) and z = 0.5 (upper curves) 
is shown as a function of the angle. The three curves at each redshift are as in figure 1. The 
two redshifts are chosen to span the range of redshifts that provide the peak contribution 
for source galaxies at Zg ^ I — 2. 
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Fig. 3. — The ratio of p{9) computed using the nonhnear spectrum to that with the hnear 
r = 0.25 CDM spectrum is shown. The normahzation is (Xg = 1, and Zg = 1. The sohd curve 
is for Vtm = 1, the dashed curve for Qm = 0.3 and the dotted curve for Qm = 0.3, Qa = 0.7. 
The results show the significant enhancement in p{6) due to nonlinear evolution for 6 < 10'. 
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Fig. 4. — The logarithmic contribution from wavcnumbers k to Cpp (dashed) and a'^{l/9) 
(dot-dashed) at = 15', 2; = 1 is shown. It is compared to that of (sohd), the variance 
in density on 8/i^^Mpc scale. The dependence oi is nearly the same as that of Cpp up to 
k ~ 0.5/iMpc~^, after which it remains positive while Cpp goes negative. 
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Fig. 5. — The dependence of p on ^ is shown for three cosmological models, with = \ 
and a r = 0.25 CDM spectrum. The thick sohd curve is for = Ij the dashed curve for 
= 0.3 and the dotted curve for = 0.3, I^a = 0.7. The thin sohd curve is for the 
hnear spectrum with = 1. The three panels show three alternative normalizations of the 
power spectrum. In the case of COBE normalization fi^ = 0.4 has been used for the open 
and A— models; see text for details of the other parameters for the COBE models. 
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Fig. 6. — The square root of the dimensionless elhpticity power spectrum a{l/9) is shown 
for the same three cosmological models as in figure The left panel uses cluster-abundance 
normalized spectra, while the right panel uses COBE normalized spectra. The thin solid 
curve shows (ril/O) for the VLm = 1 model computed using the linear power spectrum. 



-34- 



0.1 



0.06 
0.04 



0.02 







1 1 1 1 1 1 1 1 1 1 1 1 1 

\_ 

^\ Cluster CTg 

\ 




\. \^ 




\ \ 

























1 10 



9 [arcmin] 



Fig. 7. — The effect of tlie sliape of tlie power spectrum on p{9) is sliown by comparing 
tfie r = 0.25 CDM spectrum (solid curves) witli tlie V = 0.5 CDM spectrum (dot-dashed 
curves). The upper set of curves is for Zg = 2, and the lower set for Zg = I- The F = 0.5 
spectrum has more small scale power and therefore predicts larger p at small 9. 
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Fig. 8. — The effect of increasing source redshift Zg on p{6) is shown for 6 = 2' (upper 3 
curves) and 6 = 15' (lower curves). The three models shown are the same as in figure ^ and 
use cluster-abundance normalized spectra. The A— models, shown by the dotted curves, 
predict the fastest growth with discussed in the text. 
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Fig. 9. — Skewness S = {k^) / {k^)^^"^ as a function of 9 for the flat model (solid curves), open 
model (dashed curves) and cosmological constant model (dotted curves). The upper set of 
curves is for Zg = 0.5, middle for Zg = 1 and bottom for Zg = 2. Over most of the regime 
of interest the nonlinear effects are important. Note that for 6 ^ 15' perturbation theory 
underestimates the skewness of the density (Colombi et al. 1996), so that on arcminutes 
scales S could be nearly a factor of 2 larger than shown in this figure. 
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Fig. 10. — The quantity {6/5')p'^{6) 10^ is shown for the three models of figure ^. This 
provides an estimate of the signal to noise in the measurement of in a one square degree 
field with 2 x 10^ galaxies, as obtained from a typical observation with a limiting magnitude 
/ = 26. Nonlinear evolution causes the curves to be much less sensitive to 6 in the range 
2' < ^ < 2° compared to the linear prediction, which shows a significant peak around 6 = 1°. 
The linear curves for the three models are shown by the thin curves. 
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Fig. 11. — Photon propagation relative to the source-observer hne: a photon is emitted at the 
source and observed at the observer's position in the direction relative to the unperturbed 
source-observer direction (which is curved on an Euchdean plane because of background 
curvature). A deflection at x' by 5a — — 2V_l05x leads to the transverse excursion at x 
given by 5x_l(x) = r(x - x!)5a. 



